Influence of correlations on transitive electron-phonon couplings in cuprate 

superconductors 
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We investigate a model for the Cu02 plane of high-Tc superconductors where the charge carriers 
are coupled to Aig and Big symmetric out-of plane vibrations of the oxygen atoms in the presence of 
local Hubbard correlations. The coupling is implemented via a modulation of the hopping integral 
and we calculate the renormalization of vertex and pairing scattering functions based on the time- 
dependent Gutziller approximation. Contrary to local electron-phonon couplings we find that the 
transitive coupling can even be enhanced by correlations for certain momenta and symmetries of the 
vibrations. While this effect may be important for certain properties, we find that, with regard to 
superconductivity, electron-electron correlations still generically lead to a suppression of the pairing 
correlations. Our results allow for an estimate of correlation effects on the electron-phonon induced 
pair scattering from weak electron-electron interactions up to the Mott regime. For onsite repulsions 
relevant to cuprate superconductors our calculations reveal a significant contribution of Big phonons 
to d-wave superconductivity. 

PACS numbers: 71.27.+a, 74.20.-z, 74.25. Kc, 74.72.-h 



I. INTRODUCTION 

The search for the new superconducting cuprates by 
Bednorz and Miiller^ was inspired by the idea that in sys- 
tems containing Jahn- Teller (JT) ions the strong interac- 
tion of electrons with local distortions may lead to high 
transition temperatures. While meanwhile the contribu- 
tion of a JT mechanism to superconductivity is controver- 
sial, the coupling of electrons to the lattice in cuprates 
is evidenced in numerous experiments. These include 
among others the observation of an isotope effect in scan- 
ning tunneling microscopy (STM) data^ and in Tc away 
from optimal doping;^ the softening of certain phonons 
under doping)^"— and the existence of a critical CUO4 
octahedra tilt angle for the existence of superconductiv- 
ity in lanthanum cuprateSi^ Also the kink in the elec- 
tronic dispersion as seen in angle-resolved photoemission 
spectroscopy (ARPES) experiments has been attributed 
to strong electron-phonon (el-ph) interactions^ which in 
Bi2Sr2CaCu208+5 compounds is supported by the ob- 
servation of an isotope effect in the corresponding energy 
scaleiS (cf. also Ref. [ll| and references therein). How- 
ever, density functional theory (DFT) usually yields an 
electron-phonon coupling constant which is too small in 
order to explain the high transition temperatures of the 
cuprate materials. This becomes even worse when the 
effect of electronic correlations is considered which are 
underestimated within the DFT method. Basically, cor- 
relations suppress the charge fluctuations and thus the 
scattering of electrons from phonons, which leads to a 
further reduction of the el-ph interaction (for a review 
cf. Refs. [Hill. 

In terms of Hubbard type models the el-ph interaction 
is usually investigated for couplings which arise from an 
expansion of on-site energies in terms of some vibrational 



coordinateS'i^ii^"— so that in general both interactions, 
electronic and electron-phonon, are local or at least have 
a dominant local contribution. 

The interplay of el-el and el-ph interactions leads to a 
screening of the latter over the whole Brillouin zone which 
is most pronounced at large momentaj ^^i^^ For sizeable 
correlations the maximum in the phonon self-energy is 
then shifted from the Peierls wave- vector to q = 0, which 
for strong electron lattice couplings eventually induces a 
phase separation instability. 

Within a slave-boson technique the influence of corre- 
lations onto the el-ph interaction has also been investi- 
gated within a three-band model where the phonons have 
been coupled locally to the copper and oxygen densityj^ 
It was found that generically the renormalization of the 
el-ph vertex depends on vpq/uj {vp is the Fermi veloc- 
ity) . It turns out to be strongly suppressed if this ratio is 
^ 1, i.e. in the static limit, while it is almost unaffected 
in the dynamical limit. 

While the interplay between correlations and local 
(Holstein-type) el-ph interaction is quite well understood, 
the situation for transitive couplings is less clear. In a lat- 
tice model this type of coupling arises from the modula- 
tion of hopping integrals due to lattice vibrations. Vertex 
renormalization effects have been studied in the context 
of a tJ-model in Ref. [10 and for the J7 — >■ 00 Hubbard 
model within a slave-boson scheme in Refs. [21,22. We 
are not aware of any similar investigation for transitive 
couplings within an approach, which can cover the pa- 
rameter space from weak coupling up to the Mott regime, 
which is the task we want to accomplish in the present 
paper. 

In cuprates, transitive electron-phonon couplings are 
believed to be of significance for the so-called breathing 
and buckling modes which correspond to in- and out-of 
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plane motions of the oxygen atoms and which contribute 
most to the pairing interaction arising from phononsj^i^ 
For example, the inclusion of a screened Coulomb in- 
teraction within the LDA-U method^^ greatly enhances 
the el-ph interaction of the half-breathing mode whereas 
the influence on the full-breathing mode is much weaker. 
An extensive analysis of the strength and material de- 
pendencies of theses phonon modes, including also apical 
oxygen modes, has recently been performed in Ref. [26l . 
It was found that the planar oxygen Big c-axis modes 
lead to the strongest enhancement of d-wave supercon- 
ductivity. Moreover, upon including long-range Coulomb 
interactions, which in a layered system lead to low energy 
plasmon excitations fipi, the coupling can even be over- 
screened, i.e. the Big vertex increases with increasing 
^Ipi . Again the effect of strong local correlations on the 
coupling of charge carriers to Big modes has not been 
considered in Ref. [2^ and is the purpose of the present 
paper. 

Within a one-dimensional model we have previously 
studied the influence of a transitive el-ph interaction 
on the phonon self-energy Y.phiq, uj) .^"^ Most interest- 
ingly it was found that I]p;i(q, w) can be enhanced for 
small momenta upon switching on the Hubbard inter- 
action whereas for large momentum Sp?i(q, w) becomes 
suppressed with U similar to the case of a local (Holstein) 
el-ph interaction. 
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FIG. 1: (Color online) Geometry of the buckled Cu02 plane 
with Cu (small, red) and oxygen (large, blue) ions. Electronic 
degrees of freedom are on the Cu sites and vibrational degrees 
on the oxygen sites. Notice that the configuration where oxy- 
gens on the x-bonds are distorted in the negative z-direction 
leads to the same model hamiltonian as Eq. IT}. 



(Bi2212) whereas the static Q = (tTjTt) deformation has 
its realization in the tilting of the CUO4 octahedra in the 
low-temperature orthorhombic phase of La2-a;Sra;Cu04 
(LSCO). Fig.[T]shows the direction of positive oxygen dis- 
placements Ui on top of the static deformation. The 
model which we introduce in the nect section is valid for 
both cases, Q = (0,0) and Q = (vr,7r). However, for the 
latter distortion one has to take into account a staggered 
positive direction of Uix(^yy Therefore our investigations 
hold for a broad class of cuprate materials. 

For the sake of clarity it should be stated that the inter- 
play between correlations and el-ph interactions strongly 
depends on the model under consideration. For exam- 
ple, a coupling to local density fluctuations in the three- 
band model for cuprate superconductors may transform 
into a transitive coupling when one derives an effective 
low-energy one-band model by eliminating the high en- 
ergy degrees of freedom (cf. e.g. Ref. [2y). In our pa- 
per we are concerned with a one-band hamiltonian with 
transitive electron-phonon interactions which can be con- 
trasted with analogous calculations for a local Holstein- 
type coupling in Ref. [l^ Nevertheless one should keep 
in mind that our results may also mimic the local el-ph 
coupling in related multiband models. 

Our paper is organized as follows. In Sec. |TT] we out- 
line the formalism of our calculation which is based on 
the time-dependent Gutzwiller approximation. Sec. IIIII 
presents the results of the renormalization of the transi- 
tive EPI due to correlation effects and we conclude our 
investigations in Sec. IIVI In both Sec. Ulland Sec. IIIII we 
additionally discuss a one-dimensional model for illustra- 
tive purposes. 



II. FORMALISM 

We use an effective one-band hamiltonian where the 
hopping between neighboring sites is modulated by the 
displacement of adjacent oxygen atoms from their dis- 
torted equilibrium positioui^ The geometry of the dis- 
torted Cu02 plane is sketched in Fig. [T] Similar models 
have been investigated in Refs. i3I,.32, . 

The hamiltonian reads 



Motivated by this finding we investigate in the present 
paper the interplay between electronic correlations and 
a transitive el-ph interaction with regard to the Cooper 
scattering amplitude. We exemplify this interplay for the 
buckling mode where electrons couple to out-of-plane (z- 
direction) motions of the oxygen atoms. For a planar 
Cu02 system, due to symmetry, this kind of coupling 
would be quadratic in the z-coordinates;2§, However, for 
vibrations around static ZQ-displacements of the oxygen 
atoms one recovers a linear el-ph interaction. 

The static configuration of a frozen-in Q = (0, 0) mode 
where all oxygen atoms are displaced along the same 
direction (cf. Fig. [T|) corresponds to the buckling of 
the plane in YBasCuaOr (YBCO) and BiaSraCaCuzOg 



tr ^ Ui,a{c\^Ci+a,a + h.C.) + ^ 
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Here c^J creates (annihilates) an electron on a square 
lattice and the vibrational degrees of freedom are sitting 
on the bonds. The parameter tij is the hopping ampli- 
tude between sites i and j. U denotes the local Hubbard 
repulsion and Ui^a with a — x,y is the displacement from 
equilibrium of the (oxygen) atom at position i + a/2 (cf. 
Fig. [T]). gtr is the coupling parameter between lattice dis- 
placement and electronic nearest-neighbor hopping. The 
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last term in Eq. ([Ij corresponds to kinetic and elastic 
energy of the oxygens. 

Our approach is based on the time-dependent 
Gutzwiller approximation (TDGA)'^* which gives an ac- 
curate description of charge fluctuations of the Hubbard 
model^i^ and which previously has been succesfully 
applied to the study of vertex renormalizations in the 
Hubbard-Holstein model.^* The idea is to evaluate the 
energy of Eq. in the presence of an arbitrary lattice 
distortion and charge fluctuation within the Gutzwiller 
approximation (GA)^"— and to expand the functional 
up to second order in the lattice and electronic den- 
sity fluctuations. 

The Gutzwiller variational wave function has the form 
I^Pg) = Pg\SD) where Pg is the Gutzwiller projector and 
\SD) a Slater determinant. For \SD) we use a state with 
an arbitrary charge profile. We define the associated one- 
body density as pija = {SD\c^j^Cia\SD) . Since magnetic 
and charge excitations decouple we can restrict to Slater 
determinants with pij^ = pij^. 

Our expansion below will be in terms of the density 
fiuctuations of the SD which physically describes quasi- 
particles and not physical electrons 4^ The two densities 
coincide in the onsite case pua while for off-diagonal den- 
sity the real particle densities are given by Zia-Zja-Pija 
where Zia- denote the Gutzwiller renormalization factors. 

The quasiparticle density fluctuations that enter into 
the problem are 



5X„ = 



^ V(k, q)6pk+q,k,a 

k,(T 

I 1 

eO(k) + e"(k + q) 

2 C0S(fc2; 



E 



\ 2cos{ky + ^) 



'k-\-q,k,(7 



(2) 



with e'^(k) being the Fourier transformed of the hopping 
tij. In the vector Eq. ^ the first two components SX^ 
and SX^ correspond to local and transitive density fiuc- 
tuations, the latter being generated in the expansion of 
the kinetic energy. The SX^ and 6X^ components are 
also transitive fluctuations but arise from the expansion 
of the el-ph coupling. 
The expansion results in 
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ph ' 



(3) 



where the first term describes the interaction between 
electronic density fiuctuations 



^(2) = _Ly 



5X^ 
5X1 



and the elemens of the interaction kernel 



Bq Cq 



(4) 



(5) 



are defined in the appendix. 



The interaction between electronic density and lattice 
fiuctuations — l/"\//VX]i ^^P [~**3^i+a/2]'"i,Q 
scribed by the second term 
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2gtrKBq cos(y) 

2gtrKCqC0s[^-) 
2 A 
2A 



(6) 



(7) 

(8) 
(9) 
(10) 



and K = j^'Y^i^ ^cos.{kx){''^k,a) ■ zq denotes the hop- 
ping renormalization factor in the GA which generi- 
cally decreases with U and at half-filling vanishes at the 
Brinkmann-Rice transition Ubr-^'^ 

In the derivation of Eqs. (j4l)-([6]) [and the fol- 
lowing Eq. (ITTt ] we have already 'antiadiabatically' 
eliminatedi^i21i2ii^ the double occupancy fluctuations 
which arise in the expansion of the GA functional within 
the TDGA prescription. For the present model the corre- 
sponding formal procedure is analogous to that described 
in Ref. 27. 

Notice that on the GA level the coupling to the lattice 
fiuctuations is reduced ~ Zq by the correlations. How- 
ever, the TDGA induces additional couplings to local 
{6X^) and transitive (SX^) density fiuctuations and we 
will show below that in certain cases this will even lead 
to an enhancement of the electron-phonon vertex. 

Finally we obtain for the lattice energy 



(11) 



9, a/3 



and it turns out that the elimination of the double occu- 
pancy introduces a novel renormalization and a mixing 
of the X, y-buckling modes: 



[ilfr = colS^,~^'qCos{^)cos{^) (12) 
1 



2 



— {2gtrrK^KpCq. 



"q ~ i\,T'^"i''^'^ ) ^^a^^p^q- (13) 
The eigenmodes are obtained from the transformation 



yBig 



(14) 



with 
U 



COs2(f )+C0s2(f ) 



cos(^) cos(^) 
-cos(^) cos(4) 



and correspond to out- (i.e. Big symmetry) and in- 
(i.e. Aig symmetry) phase oscillations of the x/y oxy- 
gen atoms. 
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The eigenfrequencies read as 



q 



+ C0S2(|)). 



Thus the frequency of the A^g mode softens due to 
the ehmination of the double occupancy fluctuations as 
a natural consequence of the interaction between both 
bosonic degrees of freedom. There will be an additional 
softening due to phonon self-energy effects. In the fol- 
lowing we plot results in dependence of the coupling pa- 
rameter A = gl^/ujQ using the hare phonon frequencies. 
One should therefore keep in mind that the inclusion of 
the correlation induced phonon softening would give an 
additional (small) enhancement in the Aig channel. 

In order to compute vertex corrections and the effec- 
tive interaction between quasiparticles we introduce cor- 
relation functions Xq which describe the response of the 
quasiparticle (Slater determinant) density fluctuations 
[Eq. ([2])] to an external perturbation 5X\^ via 




It is convenient to represent the non-interacting sus- 
ceptibilities as 



where e.g. the 2x2 susceptibility matrix x\^'^ corre- 
sponds to the submatrix of Xij{^) for iii = 1; 2 and the 
definition of the other block matrices follows analogously. 
Notice that here and in the following we denote 4x4 
(2 X 2) matrices with a double (single) underscore. 

The lattice induced (bare) interaction is obtained by 
assuming that the phonon frequencies are larger than 
the frequencies of the electronic excitations of inter- 
est. Therefore we can eliminate the lattice eigenmodes 
for a frozen electronic configuration, dE'^^ / dV^ = 
(s = Big, Aig). One obtains in addition to Eq. ^ a 



lattice induced (bare) interaction between density fluc- 
tuations 



1 

'2N 



E 



1 



q,s,al3 1 



-11 



(15) 



Using the TDGA, this interaction together with the bare 
Hubbard part Eq. Q can be dressed by the electronic 
density fluctuations similar to conventional RPA theory. 
We define the 'undressed' interaction from the sum of 
Eqs. (I5TT51) 





y'^(q)=F-(q)+( - 



(16) 



which upon resumming the RPA leads to the effective 
dressed interaction between quasiparticles 



Ty^//(q) = 1/0 (q) 



l-X,V"(q 



(17) 



Finally we can use Eq. ([Tt)) to define the Cooper pair 
scattering amplitude 



r(k,q) = ^ ^^.^■^(q)F,(k,q)F,(-k,-q) 



(18) 



and the Vi(k, q) have been defined in Eq. We will 
usually take both k and k -I- q on the Fermi surface and 
restrict to the static limit [lu — {)). 

We can expand Eq. ((T5)) up to order ~ g^^ in the 
phonon contribution. The coefficient of ^ gf^. defines 
the second order el-ph contribution to the pair scattering 
amplitude as 



r?Hk,q) 



■^^^■^^^4F,(k,q)T/,(-k,-q) 



(19) 



This scattering amplitude does not include the direct 
scattering amplitude due to el-el interactions. It includes 
only el-ph interactions which are, however, renormalized 
by the el-el interactions. 

Here the dressed vertex functions are given by 



l-Wa 



,ee,0 



Wa 



1{ 



ei-,0 



Jq 



(20) 



From Eqs. (I20I19P it turns out that the original cou- 
pling between buckling modes and density fiuctuations 
(79'^' 7g'^)' Eq. dni), does not get renormalized by the 
TDGA but it is only reduced via the factors already 
present on the GA level. This also implies that within 
standard (Hartree-Fock based, i.e. = 1) perturbation 



theory a transitive el-ph interaction remains unrenormal- 
ized up to second order in the coupling constants. On the 
other hand, the TDGA can even induce an 'antiscreening' 
of the local and transitive fiuctuations, leading to an en- 
hancement of the 7*'*' as compared to the bare couplings 
7*'" (cf. next section). 
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Illustrative one- dimensional model 

In order to better illustrate the influence of correlations 
on the el-ph interaction within our TDGA approach, we 
will consider also a simpler one-dimensional model (cf. 
Fig. [5]). The result of this and the two-dimensional sys- 
tem will be presented in Sec. IIIII 
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FIG. 2: One-dimensional CuO model where the oxygens (light 
shaded circles) oscillate perpendicular to the chain. The hop- 
ping between nearest-neighbor Cu (full circles) atoms i and 
j -f 1 is then linearly expanded in m. 

The transformation of the previously developed for- 
malism to the one-dimensional model is straightforward. 
Since there is only one vibrational degree of freedom the 
relevant fluctuations are 

k,cr 

= J2 I + + 9) I Spk+,,k,. (21) 

k,a \ 2cos{k + q/2) I 



with a dispersion 
one obtains 



-2tcos(fe). For the bare vertices 



= 2gtrKBgCos{q/2) 
7^ = 2gtrKCgCos{q/2) 



(22) 

(23) 
(24) 



with K — —sinikp)- 



III. RESULTS 
A. One-dimensional case 

Before turning to the two-dimensional cuprate model 
it is instructive to analyze the model for a linear CuO 
chain, where the oxygens oscillate perpendicular to the 
chain direction. Fig. |3] shows the elements of the elec- 
tronic interaction kernel Eq. ([5]) for q = Q. Aq^Q cor- 
responds to the interaction between local density fluc- 
tuations. Close to half- filling it develops a pronounced 
maximum at the Brinkmann-Rice transition Ubr. where 
it actually diverges at exactly n ~ \. The coefhcient 
Bq^Q mixes local and transitive density fluctuations and 
vanishes at half- filling. However, away from n = 1, \Bq\ 



increases upon approaching U br and also stays finite be- 
yond this value [cf. Fig. [31[b)]. Finally, the coefficient Cq 
[Fig. Wis)] corresponds to the interaction between tran- 
sitive density fluctuations. It also develops a maximum 
close to Ubr^ but vanishes in the J7 — >■ cx) limit. 




FIG. 3: (Color online) (Color online) Elements of the elec- 
tronic interaction kernel Eq. ([5| for g = as a function of U 
and different fillings in the one-dimensional CuO model. The 
vertical line at U/t — 32/tt corresponds to the Brinkmann- 
Rice transition. 

Let us now turn to the discussion of the renormal- 
ized vertex functions via the scattering amplitude F^^-'/zq 
which is shown in Fig. 01 Notice that we divide by the 
GA renormalization factor which originates from the den- 
sity of states renormalization N* = Nq/zq within the 
GA. Indeed, upon plotting F^^-'/zq we take into account 
the interaction in both quantities. We point out that the 
bare coupling on the GA level is ~ gtrZ^, which for the 
'bare' scattering amplitude would result in a reduction of 
F^/zq ~ g'ir'^Q upon increasing U/t. 

In the one-dimensional model we have just two Fermi 
points so that the scattering of pairs on the 'Fermi sur- 
face' is restricted to g = and q — 2kF, respectively. 
Consider first the q — case which is shown in Fig. 
Ufa). At half-filling the scattering amplitude is identi- 
cally zero since in this case the chemical potential is at 
Ep = so that the only non-vanishing fluctuations are 
the local ones [cf. Eq. ([21]) where for q -> the fluctua- 
tions are restricted to k — hp = ±7r/2 and e'^{kp) = 0]. 
Although these are coupled to the vibrations via the co- 
efficient Bq [cf. Eq. ([22l) ]. this latter also vanishes at 
n — 1, leading to the vanishing of the scattering ampli- 
tude. Slightly below half-filling (dashed line in Fig. 0^.) 
one obtains an increase in F'^-*/zq which can be shown 
to be due to the behavior of the renormalized coupling 
lq=o ~ Bq=o/{l + Aq=oN*) to the transitive density fiuc- 
tuations. It is the concomitant increase of \Ba 
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FIG. 4: (Color online) Second order pair scattering amplitude 
-r(2'(k,q)/(4A2g) for the one-dimensional model as a func- 
tion of U. a) Scattering with q = at k = —hp; b) Scattering 
between ik^ with q = 2kF. The vertical line at U/t = 10.19 
marks the Brinkmann-Rice transition. 



Ubr together with the reduced screening from local den- 
sity fluctuations (~ 1/(1 -I- Aq^oN*) for U/t > Ubr, 
which leads to the peculiar behavior of F^^^/zq slightly 
below half- filling and q = 0. Remarkably, close to half- 
filling the scattering amplitude is amplified by large in- 
teractions and it becomes many times larger than the 
scattering amplitude in the non interacting case. This is 
just the opposite of what happens in the case of Holstein 
like couplings. 

For larger values of the doping the bare coupling j^^q 
dominates the value of F'-^-'/zq for small values of U/t 
which therefore acquires a minimum between U/t — 
and UBR./t- 

In the case oi q = 2kF scattering the behavior is sim- 
pler. At half-filling {kp = 7r/2) it can be seen from Eqs. 
([221 1221) that only the bare coupling (gtr^o) survives. For 
n = 1 this is not affected by the screening from density 
fluctuations. The n = 1 result in Fig. 2)3 (solid line) 
therefore simply reflects the behavior of Zq as a func- 
tion of U/t and vanishes at Ubr ~ 32t/7r. Away from 
half-filling the additional screening leads to a further re- 
duction of the scattering amplitude which continuously 
decreases with increasing U/t even beyond Ubr- 



B. Two-dimensional model 

From the analysis of the ID model we have learned that 
within our TDGA analysis, correlations can enhance the 
vertex for small momentum transfer. On the other hand, 
correlations generically reduce the coupling for large q 
scattering. We now investigate whether this result per- 
sists for the model in two dimensions, introduced at the 
beginning of the previous section. 

In the following figures we discuss the renormalized 
vertex function from the second order pair scattering 
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e/ji 

FIG. 5: (Color online) Scattering function — r°(kF, q ~ 0)/zq 
(in units of 8A) for various values of U /t and zero momentum 
transfer of pairs on the Fermi surface. Coupling to Aig (a) 
and Big (b) modes. Parameters: S = 0.1, t' /t = —0.4. 



amplitude T^'^\kF,(lF)/z^ [cf. Eq. ^] in units of 
— 8A which makes this quantity dimensionless and pos- 
itive. The factor of 'eight' originates from the defi- 
nition of the bare coupling Eqs. ([51 fH)) which (for 
q = 0) is given by t/tr^oV^ [cos(fcx) — cos(A;3^)] for Big- 
and [cos{kx) + cos(fcj,)] for Aig symmetry so that 

%g1j./ujQ = 8A corresponds to the maximum scattering 
amplitude for U/t — 0. 

In Fig. [S]we show the second order scattering function 
Eq. (jl9p for pairs on the Fermi surface and zero momen- 
tum transfer, i.e. T^^\]^F^^ = 0). For Aig-symmetry 
the coupling function is ^ cos{kx) + cos{ky) so that the 
U — scattering is largest around the nodes. For a 
Big-type coupling ~ cos{kj;) — cos(fcy) the U — scat- 
tering is maximized around the antinodes. This inter- 
action becomes continuously suppressed upon increasing 
U/t. Remarkably the vertex function from the coupling 
to Aig modes is instead enhanced around the antinodal 
regions while it is also suppressed around the nodes. As 
in the one-dimensional case, the enhancement by the in- 
teraction contrasts the corresponding findings for local 
couplings. 

Figs. [HI [71 [5] display F(^)(k, q)/2o for a fixed momen- 
tum k + q of the scattered pair on the FS (cf. inset). 

For the scattering from the Aig-type mode the generic 
trend is a correlation induced enhancement of the ver- 
tex function when the pairs are scattered between the 
antinodal regions. In all other cases correlations gen- 
erally suppress the scattering from the Aig mode. Since 
the enhancement occurs only for those regions in momen- 
tum space k w ^antinodes where the vertex is small, the 
overall effect is thus a suppression of the coupling upon 
increasing U/t. 
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FIG. 6: (Color online) Scattering function — r*-^'(kF, q) for 
pairs on the Fermi surface as indicated in the inset. Cou- 
pling to Aig (a) and Big (b) modes. Dressed (solid) and 
bare (dashed) vertex along selected symmetry directions in 
the Brillouin zone. S = 0.1, t' /t = -0.4. 



In the case of scattering from the Big mode, one ob- 
serves a suppression for ah momenta k and k+q. In fact, 
as in the one-dimensional case a correlation induced en- 
hancement of the vertex is due to the coupling to transi- 
tive density fluctuations (i.e. SX^). It turns out that this 
coupling vanishes for Big symmetry leading to a generic 
suppression of the coupling in contrast to the Aig mode. 
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FIG. 7: (Color online) Scattering function — r''(kF,q) for 
pairs on the Fermi surface as indicated in the inset. Cou- 
pling to Aij (a) and Big (b) modes. Dressed (solid) and 
bare (dashed) vertex along selected symmetry directions in 
the Brillouin zone. S = 0.1, t' /t = -0.4. 
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FIG. 8: (Color online) Scattering function — r''(kF,q) for 
pairs on the Fermi surface as indicated in the inset. Cou- 
pling to Aig (a) and Big (b) modes. Dressed (solid) and 
bare (dashed) vertex along selected symmetry directions in 
the Brillouin zone. S = 0.1, t' /t = -0.4. 

In order to quantify the influence of correlations on su- 
perconductivity we calculate the average of the scattering 
amplitude over the Fermi surface, weighted by the sym- 
metry function for different gap symmetries. We consider 
the following function 



(r) 



N 



N 



(25) 

and 7k specifles isotropic s-wave (is), anisotropic s- 
wave (as) and d-wave (d) symmetry: 



7k 



(cos(fc^) 
(cos(fc^) 



(is) 

cos(fcj,))/2 (as) 
cos(fc,))/2 (d). 



For the scattering amplitude r(k, p— k) we flrst use the 
second order expression Eq. (fT9)) in order to study the 
influence of vertex corrections on SC separately (cf. Fig. 
ini). Then the fully dressed expression Eq. ([T^ is used 
which, besides the repulsive Hubbard interaction Eq. 
contains all higher order contributions to the screening 
beyond Oigl) (cf. Fig. [TO]). 

Fig. [S] displays the average second order el-ph scatter- 
ing amplitude for different dopings. Notice that a nega- 
tive sign of (F) signals an attractive interaction whereas a 
positive value signals the suppression of the correspond- 
ing SC order (Ref. Houses the opposite convention). The 
largest coupling is obtained for Big symmetry where the 
isotropic s-wave pair scattering is more attractive than 
that for d-wave pairs since we are analyzing only the 
renormalized el-ph interaction without the direct effect 
of the Coulomb repulsion on the scattering amplitude. In 
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FIG. 9: (Color online) Average of the second order scattering 
function (r?'(kF,q)) (cf. Eq. ^) over the Fermi surface 
weighted with d-wave, anisotropic (as) and isotropic (is) s- 
wave symmetry functions. 5 = 0.01 (upper panel), 5 = Q.l 
central panel, 5 = 0.2 lower panel; t' /t = —0.4. The vertical 
bar in the panels marks the Brinkmann-Rice transition for 
the half-filled system. 
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FIG. 10: (Color online) Average of the full scattering function 
(rs(kF,q)) (cf. Eq. (|25p ') over the Fermi surface weighted 
with d-wave, anisotropic (as) and isotropic (is) s-wave sym- 
metry functions. 5 — 0.01 (a), 5 = 0.1 (b), and S — 0.2 
(c). t' /t = —0.4, A = g'^/i^o = 0.5. The vertical line in the 
panels marks the Brinkmann-Rice transition for the half-filled 
system. 



both channels, Big and Aig, and for all symmetries the 
average scattering from the modes continuously decreases 
with increasing U/t due to the reduction of the respective 
vertices. In the parameter range U br it turns out 
that (r^^)) is not much affected by doping. On the other 
hand one finds that for U jt beyond the Brinkmann-Rice 
limit and doping close to half-filling (Fig. [3^) the scat- 
tering amplitudes becomes negligibly small but start to 
increase for finite doping. 



The average scattering from the Aig symmetric mode 
is small for both s- and d-vifave channels since (F) is 
only large when the form factor of the SC pairs is small. 
Therefore as anticipated before the pairs cannot take ad- 
vantage from the correlation induced enhancement of the 
vertex in this case. For d-wave symmetry the effective 
el-el interaction even becomes slightly repulsive. Notice 
that in the Aig channel the anisotropic s-wave scatter- 
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ing is slightly larger than that for isotropic s-wave pairs. 
However, upon analyzing Eq. (|25l) we find that this ef- 
fect is due to the normalization (i.e. the denominator), 
whereas the numerator is larger for isotropic s-wave scat- 
tering as it should. 

Finally we show in Fig. [10] the result for the full 
averaged pair scattering amplitude for both d- and s- 
wave symmetries. The inclusion of the Hubbard derived 
quasiparticle repulsion Eq. in the interaction part 
leads to further separation between d- and s-wave sym- 
metries. Since this interaction is almost structureless in 
momentum space it mainly leads to a suppression of s- 
wave pair scattering (mostly pronounced for isotropic s- 
wave) whereas the d-wave pairs are only weakly affected. 
The latter feel the local Coulomb interaction mainly in 
the Mott-regime where at half-filling the corresponding 
scattering now also becomes repulsive. The inclusion of 
screening processes beyond 0{gfj.) via the RPA resum- 
mation starts to play a role at higher doping where it 
reduces the strong dependence of (F) on U/t. According 
to our investigations the attractive scattering of d-wave 
pairs in the Big channel dominates for U/t '> 2.5 over the 
corresponding isotropic s-wave pairing and for U/t > 5.5 
remains the only attractive pairing symmetry within our 
model. 

It should be mentioned that the GA correctly de- 
scribes band narrowing and the suppression of charge 
fiuctuations but it fails to describe the emergence of the 
magnetic J scale in the metal as the Mott insulator is 
approached from the Fermi liquid side. Therefore our 
present approach does not capture pairing due to the ex- 
change of magnetic fluctuations. Accordingly, if we ne- 
glect the el-ph interaction the scattering amplitude be- 
comes repulsive in all the channels considered. 



IV. CONCLUSION 

The aim of the present paper was to investigate the 
influence of an on-site electronic repulsion onto a transi- 
tive electron-phonon coupling within a single-band Hub- 
bard model. In contrast to similar calculations for a local 
Holstein-type couplingi^ where the electron-phonon ver- 
tex is suppressed for all momenta (but most strongly for 
large q), we have seen that a transitive el-ph interac- 
tions leads to a much more interesting interplay between 
electronic correlations and phonons. Our flnding that 
there can even be a correlation-induced enhancement of 
the pair scattering amplitude for certain momenta is es- 
pecially important with regard to the anisotropy of the 
el-ph interaction in cuprate superconductors^^. In this 
regard our approach extends to the correlated regime, 
the analysis of Refs. .2&.29 . which where based on a (elec- 
tronically) non-interacting model. 

The correlation-induced enhancement of the electron 
phonon interaction it is not followed by a similar enhance- 
ment of the overall influence of phonons on superconduc- 
tivity. Even in the present case the impact of phonons 



on superconductivity is suppressed by the Coulomb in- 
teractions. This is because the enhancement happens 
in a channel that is not beneficial for d-wave supercon- 
ductivity, namely Aig, and also because it happens in 
a restricted momentum range. Notwithstanding the ef- 
fect on superconductivity, there may be a larger impact 
on other properties which are more sensitive to the lo- 
cal, rather than the global, momentum dependence like 
the appearance of kinks on the electron dispersion due to 
self-energy effects. Also the enhancement of the scatter- 
ing function at small momentum can play an important 
role in phase separation instabilities4i 

We have presented our results as a function oi U. A 
crucial question concerns the strength of electronic cor- 
relations in the high-Tc materials. In fact, the early clas- 
sification of these compounds as doped Mott insulators 
has recently been questioned in Ref. [sl. This point of 
view agrees with our own estimate for U/t ^ 8 as derived 
from a fitting to the spin- wave dispersion of undoped lan- 
thanum cuprate o^^'^^ and which is almost half the value 
of the Brinkman-Rice transition. 

For U/t — 8 we can deduce from Fig. [TOl that the scat- 
tering amplitude for d-wave scattering from the coupling 
to Big modes is reduced by ^ 60% for (5 = 0, ^ 50% for 
6 — 0.1, and 40% for 6 — 0.2 which should be compared 
with the 20% increase when the plasmon frequency is in- 
creased to V,pi = 500meV^ Thus even the inclusion of 
(moderate) correlation effects leaves a significant contri- 
bution of phonons to d-wave superconductivity. 

In case of the Hubbard-Holstein modeli^ the electron- 
phonon vertex is determined by the local charge response 
function and it turned out that the TDGA computation 
of this quantitiy is in rather good agreement with ex- 
act diagonalization and Quantum Monte Carlo results. 
Since for the present case of a transitive coupling de- 
tailed computations from more sophisticated approaches 
are not yet available a detailed comparison concerning 
the performance of the TDGA is not possible. However, 
at least our previous investigations of charge excitations 
on finite Hubbard clusters^i^ also including the optical 
conductivity as an 'intersite' response function, suggest 
that the TDGA is a reliable approach also for the issue 
investigated in the present paper. 

Our present study may also have relevance for an un- 
derstanding of the iron-based superconductors LiFeAs 
and NaFeAs which, contrary to the other pnictide com- 
pounds, don't show magnetic order and have relatively 
low T,Mid2. which suggest that an electron-phonon mech- 
anism may play an important role. On the other hand, 
density functional theory results'*^ in coupling param- 
eters which are far too small in order to account for 
the observed transition temperature. Recently, angle- 
resolved photoemission studies^ have revealed both, sig- 
natures of optical phonons and a significant electron- 
electron scattering rate in the spectra. One of the low 
energy phonon modes in the pnictides is an iron-arsenide 
bond-stretching phonon which affects the hopping be- 
tween adjacent iron d-states in a manner similar as the 
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one analyzed here for the buckhng modes. Thus a natural 
perspective of the present approach would be its combi- 
nation with LDA computations for multiband Hubbard 
models in order to obtain a realistic estimate of correla- 
tion effects for the various phonon modes regarding their 
relevance for the SC pairing in low Tc pnictides. 
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Appendix 

The elements of the interaction kernel Eq. ([S]) read as 

= - ^ (26) 
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Notice that in Ref. the element Bq is incorrect by a 
factor of 'two'. Eq. (I?7l) reports the correct expression. 
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